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Abstract 

In the paper, the estimator for the spectral measure of multivariate 
stable distributions introduced by Davydov and co-workers are extended 
to the regularly varying distributions. The sampling method is modified 
to optimize the rate of convergence of estimator. An estimator of the total 
mass of spectral measure is proposed. The consistency and the asymptotic 
normality of estimators are proved. 
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1 Introduction 



A random IR d -valued vector X has a regularly varying distribution with char- 
acteristic exponent a > if there exists a finite measure a in the unit sphere 
S' 1 - 1 = {x | ||a;|| = l,xeR d } such that MB <= B{S d - 1 ) with a{dB) = 0, 

}^wj P {w\\ GBm>x } = a{Bl (L1) 

where L is a slowly varying function, i.e., — > 1 as a; — > oo , VA > 0. Here the 
notation || • || denotes Euclidean norm. The measure a is called spectral measure, 
and a is called simply tail index. The unit sphere S' d_1 will be simply denoted 
by S. The fact that X has a regularly varying distribution with tail index a and 
spectral measure a will be noted later by "X € RV(a,cr)". Regular variation 
condition appears frequently in the studies of the limit theorem for normalized 
sums of i.i.d. random terms, see e.g. [53] and [2], and the extreme value 
theory, see e.g. [HI]. Regular variation is necessary and sufficient conditions 
for a random Revalued vector belongs to the domain of attraction of a strictly 
a-stable distribution, if a € (0, 2), see e.g. [T]. 

There are various characterizations of the property X G RV(a,cr) (see e.g. 
|15|). We give here an equivalent definition. 

Definition 1.1. A random Mr-valued vector X E RV{a 1 a) if there exists a 
slowly varying function L such that for all r > and B € B(S) with cr(dB) = 

hm np{j^- eB,\\X\\ >rb n \ =a(B)r Q , (1.2) 

rn-oo ^ ||A || J 

where b„ — n x l a L{ri). 



It is well known that in M. d the convergence (1.2 1 are equivalent to the con- 
vergence in distribution of binomial point processes /3 n = X)fc=i^x fc /fc„ to a 
Poisson point process w a:(r whose intensity measure has a particular form |19| . 
This result is generalized to random elements in an abstract cone [S]. More- 



over, in a cone which possesses the sub-invariant norm, the convergence (1.2 1 
with a £ (0, 1) implies that X belongs to the domain of attraction of a strictly 
a-stable distribution (see Th. 4.7 [B]). 

We are interested in the problem of estimation of the tail index a and the 
spectral measure a of a regularly varying distribution. By using the relation be- 
tween the stable distributions and the point processes, Davydov and co-workers 
(see [7], [E] and [T7]) proposed a method to estimate the tail index a and the 
normalized spectral measure a of the stable distributions in K rf . The objec- 
tive of this work is to extend this method to the multivariate regularly varying 
distribution. 

Suppose that we have a sample £i» £j> • • • j £;v> taken from a regularly varying 
distribution in R d with unknown tail index a and unknown spectral measure a. 
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We divide the sample into n groups G TO ,i> ■ ■ ■ , G m ,m each group containing m 
random vectors. In practice, we choose 

n = [N r },r G (0, 1), and then m = [N/n], (1.3) 

where [a] stands for the integer part of a number a > 0. As N tends to infinity, 
we have nm ~ N. Let 

M£> =max{||e|| UeG m ,,}, i=l,...,n, (1.4) 
i the 

where the index j(m,i) is such that 



that is, denote the largest norm in the group G m ^. Let £ m .; = £j = £j( m ,i) 



11^(^)11=^- (1-5) 

We set 

M%1 = max{||e|| | £ G G^AIU,,}}, * = l,...,n, (1.6) 

(2) 

that is, M™^ denote the second largest norm in the same group. Let us denote 

M (2) - " 

_ lvl m,l q _ \ 

^raA — f-i \ i — / ^m.i 

M ( , 

m.i i=l 

and 

n - b n 

The regular variation condition ( |1.1| implies 

P{||f|| >x} = X - a cr(S)L(x) + o{x- a L{x)) as i4oo. (1.7) 

In the following we will need the stronger relation : for sufficiently large x and 
for some j3 > a 

P{||£|| >x } = C lX - a + C 2 x-> 3 + o(x- fj ). (1.8) 

Under the regular variation assumption and the second-order asymptotic 
relation (1.8), the consistency and the asymptotic normality of the estimator 
a ff were proved firstly for n — m — [s/N] in [S] and then for more general 
setting (1.3 1 in [T7]. We resume these results in the following theorem. 

Theorem A. (117$) Let . . . ,£,n be i.i.d. random M. d -valued vectors with 

a distribution satisfying (1.1) and let the numbers n and m satisfy the relation 



1.3), then 



1 „ a.s. ol 
-S n > 



n jv-s-oo 1 + a 



If the distribution of £ satisfies (1.8) with < a < f3 < oo and we choose 
n = A r2C/(i+2C)-e j m = ^ 1 /(l+20+S j 
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where £ = (/? — a) /a and e — > as N ^ oo, then as N — > oo 

•JV(0,1). 



\/n ( -S„ — ^- r 

n a + 1 




This paper focus on the estimation of the spectral measure. In [S] an esti- 
mator of normalized spectral measure &(■) = a(-)/a(S) was proposed as follows. 
We set 

m ,i = ||! m,I || , i = l,...,n (1.9) 

sm.i 



where £ m $ is defined by (1.51. Let us denote 

1 A 



(1.10) 



Random vectors # mj i, . . . , # m ,n are i.i.d. and it is proved in [8| that tx/v(0 is 
consistent considering a fixed set B. that is, VB G with a(dB) = 0, 

a N (B) ^ a(B). 

The asymptotic normality for on{B) was proved in [7]. All these relations were 
obtained under the assumption that n = m. Inspired by the work in |17| we 
modify the sampling method of regrouping and discuss the convergence rate of 



the estimator of spectral measure for the general setting (1.3 1. By finding a 



countable collection of the cr-continuity sets which is closed under the opera- 
tion of finite intersection, v 
convergence in distribution. 



tion of finite intersection, we obtain &n =$■' a as N — > oo, where =>• indicate 



Since on in (1.10 1 gives the normalized spectral measure, it remains for us 



to estimate the total mass cr(S). Note that the condition (1.1) implies 
lim x a L(x)- 1 T>{\\£\\ >x} = o(S) 

x— >-oo 

and therefore the value of cr(S) depends on the choice of slowly varying function 
L(x). Estimation of this function is discussed in |20j and |21| . Here we assume 
that the random vector £ satisfies the condition (1.1) with L(x) = 1. That 



means if a € (0, 2), the law of £ belongs to the normal domain of attraction of 
an a-stable distribution. Let us denote 



M, 



(i) 



m,2 



where is defined by (|1.4|). The proposed estimator is defined by 



a(S) 



N 



nT(l 



t > 0. 



(1.11) 



(1.12) 



4 



Example 1. We generated samples from univariate stable distribution 
with a = 1.75, cr(S) = 1 and p = ^"^g^ 1 - 1 = 0.5. The sample size is 
100, 000. We calculated the estimators p = p r as a function of r. This pro- 
cedure was repeated 50 times on the independent sets of samples. Then we 
plotted {(1 — r,p r ),Q < r < 1} where p r is the mean of 50 estimated values p r . 
Since the estimator of total mass depends on a, we present here the simulated 
results of the parameters a and p in Figure [T] The horizontal line corresponds 
to the true value of the parameter p. It seems that both plots have the opti- 
mal value of 1 — r which is around 0.4. In fact by Theorem A and the fact that 
P = 2a for a stable random variable, the asymptotically optimal value of 1 — r is 
approximately 1/3 + e. A similar result for the estimator a is given in Theorem 




Figure 1: Estimation results of parameters a (left) and p (droite) in terms of r. 
The re-axis represent 1 — r. Horizontal line represent the true value. 

Example 2. Consider the bivariate strictly stable distribution with a = 
0.75, <j(S) — 1 and the density of spectral measure defined by f(9) — || cos(20)|, 
9 £ (0, 27t). We calculated the estimators with r = 0.5. Simulated samples had 

50, 000 data vectors. The estimated parameters are a = 0.74, cr(S) — 0.99. 
The confidence intervals with level 95% are respectively (d — 0.07, d -I- 0.07) and 

(a(S)— 0.11, cr(S , )+0.13). The result of estimation of the cumulative distribution 
function (cdf ) of spectral measure is shown in Figure [2j 

The main results of the paper are contained in Section [2] and [3] The last 
section contains the proof of results presented in the previous sections. 



2 The consistency of estimators 



We assume that (,i,£,2, ■ ■ ■ ,£,n are i.i.d. random M. d - valued vectors with a regu- 
larly varying distribution. Our aim is to estimate the spectral measure a from 

the sample. The estimators <5"jv(-) and o{S) N are defined by (1.10) and (1.121. 
To establish consistency of these estimators we need two auxiliary results. 
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Figure 2: Estimates of the cdf of spectral measure in Example [2] The solid line 
is the exact cdf, the dashed line is the estimated cdf. 



Let X be a Revalued random vector, we denote G(x) = P{||X|| > x}. Let 
Yi, I2, ■ ■ • be the random variables i.i.d. with distribution function 1 — G and 
Y n _i, Y nt 2, • • • , Yn,ni Yn,i > ^n,2 > • • • > Y n _ n , the corresponding order statistics. 
Our first result extends a one-dimensional lemma in [13] (Lemma 1) to d > 1. 



It says that, for a random variable X satisfying (1.1) with <r(S) = 1, the vector 
b~ 1 (Y n i, . . . ,Y nn ,0, 0, . . .) converge in distribution to (T^ 1 ^" , T^ 1 ^", ■ ■ .). The 
multivariate version of this lemma is as follows. 

Lemma 2.1. If X E RV(a, a), then the vector b~ 1 (Y n i, ... ,Y nn , 0,0, .. .) con- 
verge in distribution in K°° to a(S) 1 ^ a (T 1 1 ^ a , T 2 1 ^ a , . . .), whereTi = ^ Xj and 

Ai,A2, ... are i.i.d. random variables with a standard exponential distribution, 
i.e. E(Xi) = 1. 

The proof is given in Section [4] 

The second result is a variant of the strong law of large numbers for a 
triangular array. 

Proposition 2.2. Let {X m i, 1 < i < n} be i.i.d. real random variables for 
each m. Suppose that the indices n and m satisfy the following relations 

n ~ N r , m ~ N x - r as N -> 00 (2.1) 

where < r < 1 is a constant and JVeN. If there exists a real number k > - 
and a constant M > such that E|X m! i| < M < 00, then 

1 ™ 

- V X mti - EX mA 0. (2.2) 



G 



Remark 1 . The convergence (2.2) holds if we replace the condition E | X m , 1 1 fe < 

oo kr 

M, k > 2 by a less restrict hypothesis J2 rrT^ 1 -^ E|X TO! i| fc < oo. 

m— 1 

The proof is given in Section [4] 

We consider at first the estimator of the total mass of spectral measure. 

Theorem 2.3 . Let £i,...,£jv be i.i.d. M. d -valued random v ecto rs such that the 
condition ( 
< t < !f, 



1.2) is satisfied with b n = n 1 /". If the condition (1.3) holds, then for 
oiS) N -a{S)-^Q, 



where o~(S) N is defined by (1.12). 

Proof. It suffices to prove the following convergence 

t ' 



1 

n Z — / 



r l 



N^oo 



(2.3) 



It follows from Lemma 



2.1 



and the assumption b n — n 1 /" that for all i and t > 
(t(S)*/ q 



as m -> oo. 



(2.4) 



Since 6„ = n 1 /", the condition (1.2 1 can be wrote as 



lim z Q P{||£|| >x} = a(S). 

Therefore there exist 8 > and a constant M > such that for x > S we have 
P{ g ^>x} = Pll/llbm 1 /", 1 /'} 

= l-(P{||e|| < m 1 /^ 1 /*})" 1 

< mP{||e|| > m 1/a x 1/t } 

< M X - a/t . 

Taking xq > S and C = xq + Mx\ °^ * we have for < t < a 



{x<x }U{x>x } 



xdP q t (x) 



< x + x P{q t m , i >x }+ V{q t ml >x}dx 

Jxq 

/■OO 

< x + Mx 1 ~ a/t + / Mx~ a/t dx 



= C. 



(2.5) 



7 



Now we choose a real number < S < 1 and 1 < 5' < 4. We set t = and 

k = &. 
such that 



Since < tk < a, it follows from (2.5 1 that there exists a constant C 



(2.6) 



In combination (2.4) (2.6 1 and Proposition 2.2 we obtain the convergence (2.3 1 



□ 



Let us consider the estimator of the normalized spectral measure defined 



by (1.101. Note that the random vectors 9 m i, ■ ■ ■ , 9 m ,n are i-i-d. in S. The 
following lemma shows the asymptotic property for each 9 mil % = 1, . . . , n. 

Lemma 2.4. Let . . . , £jy be i.i.d. M. d -valued random vectors and £ G 



RV(a,a). If 9 m ,i is defined by (1.9), then 



a as m — > oo 



for each i. 

Proof. For all Borel set B in unite sphere S such that a(dB) — 0, we have 

P{0 ro>i eB} = P{£ m ,i/||£ m ,i|| g 5} 

JTi 

= ]TP{W/IIWII 6 5,W=&} 



fc=l 



= mP{£ m /||£ m || G B,C m ,i = U} 

= mP{U/\\U\\ €B,\\U\\ > ||a||,Vfc = l,...,m-l} 

= mP{£ m /||£ m || G B, ||£ m || > b m r m -i} 

= J mP{U/Um\\ e 5, IIUII > 6ma;}P Tm _ 1 (da:). 



where P T ™_i is the distribution of r m _i = max (HCfcll^m 1 )- 

m Kfc<m— 1 



By (1.2 1 and Lemma 2.1 the last term converges to 



□ 



Therefore for each Borel set B in unit sphere S such that cr(dB) = 0, we 
have 

^B{9 m ,i) =>• ls(»7) as m ->• oo, 



where 77 is a random vector with distribution a. This yields 

El B {9 m ,i) > 



(2.7) 
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If there exists a constant r > such that n ~ N r , applying Proposition |2.2| 
for the triangular array {ls(0 m ,i), • ■ ■ , lfl(^m,n)}, we have for each fixed set 
B e B(5) with a(dB) = 0, 
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E 

i=l 



El 



B I Cm 



,l) 



N—fOD 



-> 0. 



(2.8) 



Together (2.7 1 and (2.8 1 we have the following result. 

Theorem 2.5. Le£ £, £i,...,£jv &e i.i.d. M. d -valued random vectors and £ € 
i?V A (a,(r). // <7/v(-) *s defined by (1.10) and the condition (1.3) holds, then 
VB € B(S) wito cr(&B) = 0, 



1 " 

T7 * » 



B[ym,i 



i=X 



N—too 



>a(B). 



The result is for a fixed set. A stronger convergence can be proved by an 
immediate application of the following proposition. 

Proposition 2.6. Let (S,S) be a complete separable metric space. Let {a n } be 
a sequence of random probability measures in S. If a is a probability measure on 
(S, S) such that for each set B € B(S) with a(dB) — we have the convergence 
a n (B) a -4' (j(B), then 

a.s. 

OO. 



a„ =!•' a as n 



The proof is given in Section [4] 
Corollary 2.7. Under the same assumption of Theorem\2.5\ we have 



a.s. ~ AT 

<tat =?■ a as N 



3 The asymptotic normality of estimators 



In this section we consider the asymptotic normality of the estimators a(S) N 
and (7/v- 



Theorem 3.1. Suppose that random vector^ satisfies the condition (1.8) with 
f3 > a + 1) £/ie condition (1.3) holds. If we choose 



3a - 4Q3 - 1) + - l) 2 - 8aQ3 - 1) - 7a 2 " 11 

£, it p < — a + 1, 

2a 8 



r = - — e, if /9>— a+1, 



(3.1) 



where e is an arbitrarily small positive constant, then in the following two cases 

otr 11 3 

a) < i < — Al if a + l</3< — a + 1 or B>~a + 1, 

I 4 r _ g k - 2 

3a + 2-2/3 11 „ 3 

6) < t < ^ A 1 if —a + K/3<-a + l, 

2 o 2 

(3.2) 
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we have, as N — > oo, 



i " / 1 " V 

n ^ q " hi ~ [ n ^ ^ J 
i=l V i=l / 



Af(0,l). (3.3) 



Remark 2. Fristedt, [9 , proved an asymptotic expansion for the distribu- 
tion of the norm of a strictly a-stable random vector in M. d 

G{x) = c lX - a + c 2 x- 2a + 0{x- 3a ), as x -> oo. (3.4) 

That means f3 = 2a. Therefore the condition of this theorem is satisfied if 
a > 1. If a > 8/5 the rate of convergence of estimator in C\ is close to iV 1 / 4 . 

The proof is given in Section [4] 

Before considering the asymptotic normality of the estimator of normalized 
spectral measure, we present a strong second- asymptotic relation : \/B £ B(S) 
with a(dB) = 0, 

jj^jj- € B, ||e|| > a; | = cr(.B)ar a + CaT^ + o{x~ !3 ) as x ^ oo, (3.5) 



where /3 > a > 0. Note that this condition implies the conditions (1.1 1 and 
Let us denote <t(-B) = 6, Is(^fn,i) = f) mj i,i = 1, 2, . . . ,n. Then 

1 ™ 

Z n = a N (B) - a(B) = - YWi - or(B)), 
n * — ' 

i=i 

1 n 1 n 

V"^n = /A*?m,« - cr(-B)) = — F= /A»?m,i - E?7 m ,i) + V™(E?7 m .i - 

v »=1 v i=l 

1 - 

£4 = —j= /XVm.i - E?7 m ,i), r m = Er) m ,1 - cr(i?). 
v ^=l 

We set _ 

„ a/ nZ r , 



V i=l \ i=l 



V 



1/2 ' 



then the asymptotic property of <7jv can be described as follows. 
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Theorem 3.2. Let £, £i,...,£jv & e i.i.d. M. d -valued random vectors and the 



distribution of £ satisfies the condition (3.5). If we choose 
n = AT 2 C/(i+2C)-e ) m = ^1/(1+20+^ 

where ( — min(^^, 1) and e is an arbitrarily small positive constant, then 

T N =>Af(0,i). (3.6) 

Remark 3. We can get also asymptotic normality for y/nZ n , but the vari- 
ance of the limit normal law is o~(B)(l — cr{B)) which we are estim ating . 



Remark 4. If £ is a strictly a-stable random vector in K d , by (3.4 1 we get 
(3 = 2a, thus the asymptotically optimal value of n is approximately N / 3 . The 
rate of convergence of &n{B) in L\ is close to TV 1 / 3 . 

The proof is given in Section [4] 

4 Proofs 

Preliminary remarks 

We recall the definition of a regularly varying function. We say that L is a 
regular varying function of index a at infinity (respectively at origin) and we 
denote L £ R a (respectively L € R a (0+)) if 

— > x a , as x — !• oo (x — >• 0) for all A > 0. 

L(x) 

Let X be M d -valued random vector satisfying the regular variation condition 



fL2| . We denote G{x) = V{\\X\\ > x}. Then 

nG(b n x) — > a(S)x~ a , as n — > oo, for all x > 0. (4.1) 

For positive fixed x, we choose n the smallest integer such that b n+ i > x. Then 
b n < x < b n+ i and for a non-creasing function G we have 

G(Xb n+1 ) G(Xx) G(Xb n ) 

nl , \ < n , \ < 7' fOT a11 X>0 - 

G(b n ) G[x) G(b n+ i) 



By ( |4.1| we have nG(b n ) ->• cr(S), then 

- ► ^ _a j as a; — > oo for all A > 0. 

G(x) 

We deduce that G <G R~ a , which allow us write the following equivalence 

G(x) ~ x ~ a L(x), (4.2) 

where L(x) is a slowly varying function. 

We recall a well known result on the asymptotic inverse of a regular varying 
function. 
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Theorem 4.1. (J3j/ Th. 1.5.12) Let f £ R a with a > 0, then 3g(x) £ R 1/a 
such that the following relation holds 

f(g( x )) ~ g(f( x )) ~ x as x -> oo- (4-3) 

-Here g (the asymptotic inverse o/ /J is defined uniquely up to asymptotic equiv- 
alence, and a version of g is 

f~(x) = M{y : f(y) < x}. 

We denote 

/(*) = ~ ~ x a j^r, (4.4) 
G(x) L(x) 

then f(x) £ R a with a > 0. By applying the previous theorem, we obtain the 
inverse g{x) of /(x) in the following form, 

g (x) = x 1 / a L\x) 

where the slowly varying function L" verifies the following relation 

Lixy^L^fix)) -> 1, (4.5) 

and 

L{g{x))- 1 L' i {x) a -> 1, i -> oo. 



By (4.3 1 and (4.4 1 we have 



G I g I i ) ) ~ x, x -> 0. (4.6) 



Defining the generalized inverse 

G-\x) :=M{y:G(y)<x}, 

we can prove that 

G(Cr\x)) ~x, x^O. (4.7) 

For this we choose X > 1, A > 1, S £ (0, oo), then by the theorem of Potter (Th. 
1.5.6 [3i page 25) there exists Uq such that 

A- 1 X- a - s G(v) < G{u) < AX a+5 G(v), W £ [X'^Xu], u > u Q . 

We take x small enough such that G _1 (x) > u , then by the definition of G _1 
there exists y £ [A _1 G _1 (x), G _1 (x)] such that G(y) > x, and there exists 
y' £ [G _1 (x), AG _1 (x)] such that G(y') < x. Taking G _1 (x) for u, y and if for 
v, we get 

A-iX-"- s G(y) < G{G~ 1 {x)) < AX a+s G(y'). 

Hence limsup and liminf of G(G~ 1 (x))/x are between AX a+s and A~ 1 X~ a ~ s 
as x — > oo. Taking A, X 11, we have G(G _1 (x))/x — > 1. 
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The relations (4.6 1 and (4.7 1 give immediately 



G^ 1 (x) ~ g(l/x), x^O. 
Thus we have the equivalent expression of the inverse of G(x): 
G~\x) ~ x- 1/a L\l/x) G R-i /a (0+). 



(4, 



Lemma 4.2. Let X be a R d -valued random vector, G(x) = P{\\X\\ > x}. If 
X G RV(a, a), then for each £ = 1,2,..., 

b^G- 1 ( > a(S) 1/o T~ 1/a with probability 1. (4.9) 

Proof. We recall (4.1l nG(b n x) — > a(S)x~ a which implies 

a(S) ( b n \ a 



G(x ri 



n — > oo 



where x n = b n x, x n — > oo, as n — > oo. By replacing the left term in the previous 
formula by (4.2 1, we get an equivalent expression of b n in terms of L(x): 



bn 



nL(x n ) 



X/a 



(4.10) 



Considering (4.8), we have an equivalent expression with probability 1 for 
each i, 



G- 1 



— Xj a 



v- 



r«+i 



, n — > oo 



where satisfies (4.5 1 which means 



1, n — > oo. 



(4.11) 



(4.12) 



Collecting (4.10 )-(4.12 ) we deduce that with probability 1 

XI a 



b- L G 



Xn-x 



r„ 



+i 



a{s) x/ a l igl 



n+1 



<r(S) 



l/a 



as n — >■ oo, since T n+ i/n a ' 6 ' > 1; the lemma is proved. 

n— >oo 



□ 



Proof of Lemma 2.1' It is well known that (see, e.g. [I] Section 13.6) 



Tx 



n+X 



,G~ 



r„ 



n+l 



The lemma follows from (4.9 1 and (4.131. 



(4.13) 
□ 
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Proof of Proposition 2.2' Denote Y m ^ — X m ^ — EX mi i, then the random 
variables {Y mi i, 1 < i < n\ are centered and i.i.d.. We have 

E|F m ,i| fc < E(|X TOil |+|EX m)1 |) fc < E(2 k - 1 (\X rnA \ k + \EX m . 1 \ k )) < 2 k V\X m , 1 \ k . 

It is well known (see [55]) that for k > 2 we have 



E 



<c(fc)n fc / 2 E|y TO)1 | 



where c(fc) is a positive constant depending only on k. It follows from the 
condition (2.1 1 that there exists a constant C > such that n > CN r . Hence 
for all e > we have 



1 ™ 

71 ^ 



E 



>s < 



< 



2 fc c(fc)E|X TO;1 



Co 



N^e k ' 



where c = 2 k c{k)M/Ci. Since ^ > 1, we can find a small enough positive 



number e' such that 



e' > 1. Taking e = En = N * and applying the 



Borel-Cantelli lemma, we have that with probability 1 and for N large enough 

1 



the proposition is proved. □ 
Proof of Proposition \2.6\ We denote the collection of all er-continuity sets 

by 

T) a — {B I B € B(5), o{dB) = 0}. 
Since space S is separable, there exists a countable dense set in S, denoted by 

W = {xi,x 2 , ■ . ■}, Xi e S, i = l,2, — 

We denote the open ball with centre Xi in W and radius r by 

V(xi,r) = {x \ x € S,\\x - Xi\\ < r}. 

Since for each Xi £ W the boundaries d{V(xi,r)} C {x \ \\x — Xi\\ — r} are 
disjoints for different r, at most a countable number of them can have positive 
a- measure. Therefore, there exists a sequence of positive numbers ri j as 
k — > oo for each x; such that 



h i = {V(x i ,r i k ),k = l,2,...}cV a 
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The collection L = (J Lj is countable. It is clear that for each xi € W, the 

Xi£W 

collection is a local base at point Xi for la topology S. Since W is dense in 
S, L is a base of 5. The a-algebra generated by L, denoted by cr(L), is the 
Borel-field 

Now we expand L by adding the finite intersections of members of L, we 
denote 



C = L U J p| Vi 



Vi e L, I C N, card(J) < oo 



It is clear that C is still countable and a(C) = B(S), moreover C C T) . Since 
a n (B) Q 4' a(B) for all B e B(S) and cr(aB) = 0, then W G £, 3Ay C H and 
P(Ay) = 0, such that Vw € A^ we have 



o- B (w,V)->-o-(V0. 



(4.14) 



If we denote A = (J Ay, then P(A) = 0. Moreover Vw £ we have always the 
vec 

convergence (4.14) for all V 6 £. The collection £ is closed under the operation 
of finite intersection. By Theorem 2.2 in [5] (page 14) we have a n =>■ a, Vw € A c . 
which implies a n =£' o\ □ 



3.1 



We set X„ 



Eg* 



Var(X„ 



Proof of Theorem 

and r m = ^ m - T(l - t/a)a{S) t ' a . Therefore the convergence ( |3.3| holds if we 
have the following three relations : 



n 

^E(*« 

V 2 — 1 



Mm) -A/"(0, <T 2 ) 



jV- 



P 2 
-> cr . 



(4.15) 

(4.16) 
(4.17) 

By the similar method to (2.6 1 we can prove the moments {EJf^} are 
uniformly bounded. Hence we have the following convergence 



i=l \ i=l / 



and 







- (EX„ 



a(S) 2t ^(T(l-2t/a)-(T(l-t/a)) 2 ), (4.18) 



and the Lindeberg's condition, i.e. for all e > 



Km 

n->oo no~ 



1 " 
E 



-/i m |>e,/racr m } 



= lim — 

n— >oo (J 



m,l Mm 

{\X m: i—iM m \>e^/na m } 

< lim 4^(E(X m>1 - Aim ) 4 ) 1 / 2 (P{|X m4 - Mro | > £^a m }) 1/2 



n—too <j£ 

= o. 
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The convergence ( |4.15 1 follows from the central limit theorem applied to tri- 
angular array {X m> i,l < i < n}. By Proposition 2.2 and inequality (2.5 1, if 
< t < 2p we have the following convergences 



and 



It follows 



1 " 

n i 

i=i 

1 " 

n ^ 



N- 



-> 



a.s. r , 

(j, m > o. 

JV-S-OO 



1 " (l n \ 2 

- X! ~ I ~ X! Xm ' i ) ~ ( E ^m,i _ Mm) 



N- 



->0, 



considering (4.181 we obtain (4.161. Now it remains to verify(4.17 1. 



Lemma 4.3. If £ satisfies the condition (1.8) with (3 > a + 1 and < t < I, 
then 

\r m \ < Cm-i (4.19) 
where £ = min(^, t+,3 ' a Q; ~ 1 ) and C depending only on ci, C2, a, /3 and i. 



Proof. By relation (1.8 1 if X is sufficiently large we have 

G(x) = Cl x~ Q + c 2 x~? + o(x- ), 

where c\ = u{S). It is possible to write the inverse function for small value of t, 

G _1 (t) = a(S) 1/a t-^ a + bt s + 0{t s+ ^-^' a ) 

with b = a- 1 c 2 <7{S) ( - 1 ~^l a et s = (0 - a - 1) /a. It follows that for small 5 > 
and < t < £ 

G -1 (t) - (i(S) 1/a r 1/Q = t s (6 + 0{t {f3 - a)/a )). 

We choose 5 such that \0{t^~ a ^ a )\ < \b\ (this gives us 5 = o(\b\^)), then we 
can write 

\G-\t) - a(S) 1/a t'^ a \ < 2\b\t s . (4.20) 



Let 



= {x = (xi,. . . ,x m+ i) : Xi > 0,i = 1, . .. ,m+ 1}, 
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where 5 is chosen for that (4.201 holds. By its definition and the relation (4.131 
the random variable q m ^ = M^ > i /m 1 / a have the distribution 



G~ 



m 



then 

km I = 

where 
h = 



E | G- 1 



_Ti_ 



l*/ Q - E 



t/c 



<h+h+h+h, (4.21) 



-t/a 



-t/c 



G 



i ( x i 



A" 



G 



^m+l 
-1 / x l 



exp(-S m+ i)(ix, 



'-"m+l 



tr(S) t/a 



exp(-£ m+ i)d 



h 
and 
h 



exp(-E m+ i)dx, 



Since in A we have G 1 (xi/S m+ i) < G 1 (5), for all t\ > (which will be 
chosen later) 

h < m-'/ Q (G- 1 («J)) t P j-"- 3 - > (?) < m- 4 / tt C*iE ( (4.22) 

i Tm+i J vr m +i/ 

where G\ = (G~ 1 (6)) t S~ Tl . In a similar way we estimate I4, for all T2 > 0, 

1/2 



Tm+l 



(4.23) 



where C 2 = cr(S')*/ Q (r(l - 2i/a)) 1 / 2 ^ T2 . The fu nction f(x) = x\0 < t < 1 is 
Lipschitz continuous on [0,(5). Considering (4.201 we have 



h < 



-*/<*c* 



-1 / x i 



< 2|6|C 3 TO" t/a / 

J A 

< 2\b\C 3 m- t/a E 



-1/a 



A° V ^m+1 



exp(-S m+ i)dx 



<t(S) 1/q 



r 



m+l 



exp(-£ m+ i)da 



(4.24) 



where C3 = £(min {G- 1 {5),5- 1 / a (j(S) 1 l a )) t - 1 . It is well known (see for example 
j¥|) that the m-dimensional random vector ( — 



L m+l 



m+l 



has the same 
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density as the order statistics vector of random variables uniformly distributed 
on [0,1). In particular the random variable — has the following density 

Im+l 

(suppose that m > 2) 

f m^-xY 1 - 1 , if 0< x< 1, 
^ ' \ 0, otherwise. 

Thus 

F,(-^— ) = m/ x a {l-x) m - l dx = mB(s + l,m) < AT(s + l) m - s . (4.25) 

By Cauchy Schwarz's inequality we have 

h < m- t ^ a a{S) t/a (ET- 2t/a ) 1 ^(E(T t J l a - m'/ a ) 2 ) 1/2 . (4.26) 

It remains to evaluate 

E(r^-m^) 2 = r(r»+l + 2f/a) _ t/a r( m + l + t/a) w/ 
v m ; r(m + l) r(m + l) 

Using the approximation of Gamma-function 

r(,) = ^-^(i + l + ^ + (^ 

we obtain 

E(r^ Q -m*/ a ) 2 = A m (™+l+2t/a) 2t/Q -ffi-2B m m*/ Q (m+l+Va) t/Q ^2+m 2t/Q , 

(4.27) 

where 

/ i/a \ m+1 / t/a Y 1/2 tla 
\ m + 1 J \ m + 1 J 

H 1 = l--^ + o(- i -) ct H 2 = l-—^+o(^). (4.28) 

By a simple calculation we obtain the following expansion 

* 2i 2 \ 1 (\ 

~ + — H 7 +o - 

a ar/TO + 1 \m 



^ = l-( 7 + ^| + o ( — ) , (4.29) 



Bi " = 1 l£ + 5?);nW=]. 



Considering ( 4.27 1-( 4.30 1 we have 



E(I^ Q - m'/") 2 = ^m 2 */"- 1 + aim 2 *'*- 1 ). 
or 
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Therefore it follows from (4.261 that there exists a positive constant C4 such 
that 

h < Cim- 1 ' 2 . (4.31) 



Collecting estimates (4.21 1- ( 4.25 1 , (4.31 I and choosing r\ = (ft — a — l)/a and 
r 2 = (t + fi - a - 1)J a in ( |4.22| and ( |4.23[ ), we obtain the estimate ( |4.19| with 
C = C(a, /3, c\,C2,t). The lemma is proved. □ 

Having the estimate for r m , and taking n ~ 7V r with r < y^§^ we obtain 
y/nr m — > 0. By a simple calculation we have the relations (3.1 I and (3.2). □ 



Proof of Theorem 



3.2 



We denote c 2 = b(l — b). In order to prove (3.6) it 
suffices to show the following relations, 



JV-s-oo 



>AA(0,c 2 ), 



\fnr m 0, 



(4.32) 
(4.33) 



1 n I I n 

-J2<i- -E 



f 2 
-> C . 



(4.34) 



Since < I.B(0 m ,i) < 1, the moments E|ls(6'„ 
formly bounded for all k > 0. The limit variance is 



El 



6 + r m - (b + r m ) 2 ->■ c 2 , 



if we have (|4.33|). We consider the Lindeberg's condition: for all e > 0, 

(l B (e„ M )-EI B (0 mjJ )) 2 dP 



lim 

n->oo riC; 



lim 



1 



l .,)-El B (e m ^)\>e^ic m } 

(ls(^,i)-El B (0 m a)) 2 dP 



'{|l B (e m>1 )-El B (e mil )|>£ % /Sc m } 
< lim i-(E(l B (0 m ,i) - El B (0 mil )) 4 ) 1 / 2 (P{|? ?m a - Ery ma | > e^cj) 1 / 2 



= 0. 



(4.35) 



The relation (4.321 follows from the application of central limit theorem for the 
triangular array {lB(d m ,i) ~ E1,b(6>„ m )}. 



It is easy to see that (4.341 follows from the application of Propo sition 2.2 
to triangular array (j)^ J. It remains to establish the relation (4.331. In fact, 
we have |r m | < C max.(m^ 1 , m~^~ a ^ a ) from the following lemma. Taking an 
arbitrarily small positive constant e and 



n = N 1 + 2 ^ e , m — iVi+2<+ e ; 
we get \fnr rn — ¥ 0, the theorem is proved. 



□ 
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(4.36) 



Lemma 4.4. If the condition (3.5) is satisfied, then 

\r m \ < Cmax(m- 1 ,m- (,3 ~ Q)/Q ). 
Proof. In order to prove ( |4.36 1, we need to show that 

P{9 m , i eB} = a(B) + R m , (4.37) 
with the remainder term R m = 0(max(m~ 1 , m~^~ Q )/ a )). Let us denote 



G m (x)^'P{ max ||^|| < x 

1 KKm-l 



Using the definition of 6* m it is not difficult to see that 

P{0 m ,i EB} = mJ^ pj^GB, Ui\\ > r\ G m (dr). 



Let G m (x) — G m (xm}/ a ). Assumption (3.5) implies (we suppose that cr(S) = 1) 
that for large s, 

P{M\\ > s} = s- a + Cs-P +o(s-P). 



Therefore, it is easy to get the relation 
lim G m (x) = G (x) = 



exp(— x Q ), x > 0, 
0, x < 0. 



Using (3.51 and the fact that J Q y a dGo(y) — 1, we have (4.37) with R m = 



Z)i=i R rn : i, where 

Rm,x = mf ^ B MA\ >r}dG m (r), 

R m ,2 = -a{B) f y- a dG (y), 

J O 

/CO 
y- a d(G m (y)-G (y)), 

nOO 

RmA = Cm-^/o / y-PdG m (y). 



Here s' = sm 1 / Q and we shall choose s later. It is easy to see that 

Rm,i < m(l - P{||?i|| > s}) m = m(l - fa- 1 )" 1 < me-i h , 

where h — h(m,s) = mP{||£i|| > s} > \ms~ a . We have used (jZJ) for the last 
inequality. Thus, if we choose 

l/a 



\K\nm) 
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with sufficiently large K, then we get 

R m .i = oim' 1 ). (4.38) 

Simple calculations show that 

R m ,2 = o{m- x ). (4.39) 

The main remainder term is i? mj 3 and to estimate it we must first estimate the 
difference G m (y) — Go(y). A rather simple expansion of logarithmic function 
gives the following estimates which are sufficient for our purposes. 

Lemma 4.5. ([!} Lemma 2) Let > 1 be i.i.d. random vectors satisfying 
Then for y > cmT 1 / 01 

\G m (y)~G (y)\<C(a,(3)eM~y' a )^- {p - a)/a y- p +m- 1 y- 2a ) (4.40) 

and 

sup \G rn (y) - G (y)\ = C(a, /?) m^( m -\m^ a ^ a ). 
v 

Now we can estimate the term R m ,3- Integrating by parts, we get 

\R m M=<r(B)(R% 3 +R% 3 ), (4.41) 

where 

R% 3 = s , - a \G m (s')-G ( S % 
R ( Z = a \ G ^(v) ~ G Q {y)\y-<*- l dy. 



Since s' = (if In ra) 1 ' a > Cm 1 / Q , we can use (4.401 to estimate both quanti- 
ties -R^n 3 , i — 1,2. After some simple calculations, we get 

B% 3 = o{m- 1 ), 



R% < C(a,/3)max(m- 1 ,m-^- Q '/ Q ). 
In a similar way we estimate R m ^: 

/>oo 

R mA = Cm-W- a V a y-PdG m (y) = Cm-V- a » a (RW 4 + R ( V A ), (4.42) 

J s' 

where 

/CO 
y- p dG (y), 

R^A= / y- p d(G m (y)-G (y)). 



It is easy to see that 

R&\<C(a,P) (4.43) 
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(2) 

and R m 4 can be estimated in a similar way to R m ,3'- 

R%\ < C(a,/3)max(m- 1 ,m-^- Q)/Q ). 



Collecting (4.38), (4.391, and (4.41 1-(4.44 1 we get (4.371 



(4.44) 
□ 
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